Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Spatial patterns and predictor variables vary among different types of primary producers and consumers in eelgrass (Zostera marina) beds

  • Mizuho Namba ,

    Contributed equally to this work with: Mizuho Namba, Masahiro Nakaoka

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft

    mizuhonamba77@gmail.com

    Affiliations Akkeshi Marine Station, Field Science Center for Northern Biosphere, Hokkaido University, 1 Aikappu, Akkeshi, Hokkaido, Japan, Graduate School of Environmental Science, Hokkaido University, Sapporo, Hokkaido, Japan

  • Masahiro Nakaoka

    Contributed equally to this work with: Mizuho Namba, Masahiro Nakaoka

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – review & editing

    Affiliation Akkeshi Marine Station, Field Science Center for Northern Biosphere, Hokkaido University, 1 Aikappu, Akkeshi, Hokkaido, Japan

Abstract

Eelgrass (Zostera marina) forms extensive beds in coastal and estuarine environments and provides various ecosystem functions. The aboveground part of eelgrass provides habitats for other types of primary producers such as epiphytic microalgae and for epifaunal invertebrate grazers. Because of the different sizes, generation times and resource requirements, these different types of producers and consumers may be affected by different sets of biotic/abiotic factors over multiple spatial scales. We examined the spatial variations in three functional groups of eelgrass beds (eelgrass, epiphytic microalgae and epifaunal invertebrates) and the abiotic/biotic factors responsible for these variations in three lagoons with different environmental properties at the eastern region of Hokkaido Island, Japan. The spatial scale responsible for the variation in the biomasses of the three functional groups varied, where within-lagoon variation was important for eelgrass and epifauna but among-lagoon variation was important for microalgae. The environmental predictors for the observed spatial variations also differed among the different functional groups, with variation in eelgrass biomass related to depth, nutrient and salinity, epiphytes to water temperature, eelgrass biomass and water column chlorophyll and epifauna mainly to eelgrass biomass. These results revealed that the level of importance of among- and within-lagoon environmental gradients vary in the different functional groups of the eelgrass bed community. The large-scale variation in pelagic productivity, which is tightly related to the ocean current regimes, is likely responsible for the great among-lagoon variation in microalgae. The local variations in environmental factors such as salinity and nutrients, which change with alterations in terrestrial river inputs, are likely related to the great variations in eelgrass and epifauna within the ecosystem. The observed relationship of epifauna with eelgrass biomass indicates the importance of non-trophic plant-animal interactions because epifauna utilize eelgrass as habitat. We therefore emphasize the importance of evaluating spatial variations at multiple scales to further understand the functions of coastal and estuarine ecosystems.

Introduction

Intertidal and subtidal meadows formed by submerged aquatic vegetation (SAV) are among the most common habitat types in coastal marine ecosystems [1, 2], and these ecosystems are often compared to terrestrial forests and grasslands in terms of their ecosystem functions and services [3]. Seagrass beds are one of these ecosystems and are found in shallow coastal marine and estuarine environments around the world [4]. The ecosystem functions provided by seagrass beds include the provisioning of habitats for diverse fauna and flora through the addition of physical structures to the seafloor [2, 5] and the support of high primary and secondary productivity [6], which makes them one of the most productive ecosystems in the world [7]. Primary producers in seagrass beds consist of various plant functional groups that are separated in terms of their turnover rates and light and nutrient requirements and include groups such as seagrasses, epiphytic microalgae on seagrass blades, benthic algae and phytoplankton [8, 9]. The amount of production by various primary producers differs spatially and seasonally [10, 11]. The amount of production by epiphytic and benthic algae sometimes exceeds the production by seagrass [8]. Consumers in seagrass beds are also diverse and consist of small invertebrates such as gastropods, amphipods, shrimps, and annelids, and some vertebrates such as rabbitfish, green sea turtles, manatees and waterfowl [12, 13]. The invertebrates can be categorized as epifauna or infauna depending on where they inhabit the seagrass beds, and they utilize the food sources and habitats provided by the primary producers, as most of the invertebrates are grazers and detritivores [14].

There have been many studies that have examined the variations in primary producers and consumers in seagrass beds [15, 16]. Investigations have been conducted on plant-plant interactions, such as the competition for light and nutrients between seagrass and algae [17], and plant-animal interactions, such as the control of plant abundance by herbivores [12, 18]. For epifaunal consumers such as invertebrate herbivores, the bottom-up effects from plants on animals have been investigated, and these investigations have focused on the relationships among the abundances of associated animals, microalgae and seagrass (e.g., [1921]). The investigations of the top-down effects have focused on the predation pressures on herbivores by predators such as fish [22]. Such plant-plant and plant-animal interactions are affected by the variations in multiple environmental factors such as temperature, salinity and nutrients [23, 24].

The different environmental factors that affect the ecosystem functions of coastal habitats, including seagrass beds, operate at different spatial scales. The abiotic factors of ocean currents and associated oceanographic factors (such as water temperatures and nutrient concentrations) affect the community structures of primary producers and consumers over large spatial scales; e.g., 100 km scale [25, 26]. At the meso scale (e.g., 10 km scale), the presence of temporal and spatial salinity gradients within lagoonal and estuarine ecosystems is common due to the changes in the tidal effects from the ocean and freshwater inputs from rivers [27]. At this spatial scale, the differences in land-use along the watersheds of the rivers that flow into lagoons and estuaries also cause variations in the amounts of nutrient input [28], while variations in physical properties, such as water depth, salinity and hydrodynamic conditions, impact the community structures and the amount of production [2931]. At fine spatial scales (< 1 km, for example), it is commonly observed that microhabitat heterogeneity and species interactions affect the diversity and abundance of different components of organisms in a community, as mentioned above. Different functional groups of organisms likely respond differently to these environmental variables due to differences in functional traits. Microalgae within seagrass beds have much faster turnover rates than seagrasses [32] and may have different relationships with nutrient concentrations and other environmental requirements, leading to different spatial variation patterns across multiple scales compared to seagrass. Therefore, a simultaneous investigation of different types of primary producers and consumers over multiple spatial scales is necessary to understand how these functional groups interact and structure the highly productive ecosystem. Additionally, how various environmental factors in seagrass beds affect these groups of organisms should be analyzed by comparing multiple seagrass bed ecosystems that are stretched over a large spatial scale, that are different in environmental properties at meso scales, and that environmental heterogeneity is present at fine scales for more precise evaluation of ecosystem functions.

This study aims to investigate the patterns in the spatial variation of each of the major functional groups of primary producers (eelgrass Zostera marina and epiphytic microalgae) and consumers (invertebrate epifauna) in eelgrass beds, the most widely distributed seagrass beds in the temperate northern hemisphere [15], and to determine how they are shaped by various environmental factors. First, we determine the most prominent spatial scale of the variability in the biomass of the two groups of primary producers and epifauna by comparing the among-lagoon (20–200 km; large scale) and within-lagoon (< 20 km; meso and fine scale, but hereafter we refer to this as fine scale) variation. We then explore the abiotic and biotic factors that are related to the observed variabilities in the biomasses of different functional groups. We study three lagoons on the eastern part of Hokkaido Island, Japan (hereafter eastern Hokkaido) which exhibit differences in the effect of ocean currents and the amount of river inputs. In these lagoons, eelgrass forms extensive beds and supports high primary and secondary production, which benefits the local communities through commercial fishing of seagrass-associated secondary producers such as Pandalus shrimp [3335]. Here, we hypothesized that the among-lagoon variation would be responsible for the variation in biomass if the organisms are affected more by large scale environmental factors such as the difference in water temperature, whereas within-lagoon variation would be more important if the organisms are affected by the factors associated with finer-scale nutrient or salinity gradients caused by freshwater and terrestrial inputs with river discharges.

Materials and methods

Study area

Eastern Hokkaido.

Eastern Hokkaido faces two different oceans: the Pacific on the southeastern side and the Sea of Okhotsk on the northeastern side (Fig 1). The coastal areas are characterized by the presence of multiple semi-enclosed estuaries and lagoons with saline to brackish water. The Oyashio cold current flows in the southwest direction along the Kuril Islands and Hokkaido coast [36], strongly affecting the Pacific Ocean side of eastern Hokkaido. The seawater in the Oyashio-influenced area is characterized by cold temperatures reaching < 0°C in the winter [36] and a high concentration of chlorophyll a (hereafter Chl-a) in the water column during the spring bloom from February to April ([36, 37], personal observation T. Isada). In the coastal areas, an autumn bloom from August to October has also been observed (personal communication T. Isada). The water mass of the current is rich in NO3 in the winter but is depleted in the spring when phytoplankton blooms [37].

thumbnail
Fig 1. Study sites.

(A) In eastern Hokkaido at the (B) Saroma-ko Lagoon, (C) Notoro-ko Lagoon, and (D) Akkeshi-ko Lagoon and Akkeshi Bay.

https://doi.org/10.1371/journal.pone.0201791.g001

The Soya warm current flows in the Sea of Okhotsk towards the southeast along eastern Hokkaido [38, 39]. The current dominates in May-October and is responsible for the higher salinities and temperatures than in the Oyashio current [39]. In November-April, the effects of the Soya current in the Sea of Okhotsk are weakened, and the East Sakhalin cold current becomes more influential and brings water masses with low salinities and low temperatures [40]. The timing of the spring bloom and nutrient depletion is similar to that in the Oyashio-influenced region, but the coastal areas of the Sea of Okhotsk have lower phytoplankton abundances and nutrient concentrations compared to the Pacific coastal areas [36, 37, 39]. Moreover, seasonal sea ice usually covers the Okhotsk coast from January to April [41].

For this research, no permit was required for the field sampling as the method used for the study (use of mesh bags) was exempt from the list of fishing gears that need to be declared and used with fishing permits. Moreover, the sampled plants and invertebrates did not include the commercial species or protected species.

Akkeshi.

The Akkeshi-ko Lagoon and Akkeshi Bay are located along the Pacific side of eastern Hokkaido (Table 1; Fig 1). They are connected to each other by a narrow channel that is approximately 500 m wide and 10 m deep. The surface areas of the Akkeshi-ko Lagoon and Akkeshi Bay are 32 km2 and 110 km2, respectively. The southern part of the bay is open to the Pacific Ocean. There are three major rivers (Bekanbeushi, Tokitai, and Tobai) that flow into the Akkeshi-ko Lagoon; 98.8% of the total input to the lagoon is from the Bekanbeushi River [24], and its average daily input accounts for 5.8% of the total water volume of the lagoon [42]. The watershed of the Bekanbeushi River is mostly covered by wetlands, agricultural lands, and forests. The combination of saline water from Akkeshi Bay that is influenced by the Oyashio current and the freshwater input from the rivers makes the water mass in the Akkeshi-ko Lagoon brackish and is responsible for the various physical and chemical gradients, including the presence of temporal and spatial salinity and nutrient gradients [24, 42].

thumbnail
Table 1. Environmental conditions in the summer (June to July) and fall (August to September) 2016 at the study sites in Akkeshi, Saroma and Notoro.

https://doi.org/10.1371/journal.pone.0201791.t001

The water depth of the Akkeshi-ko Lagoon ranges between 0.8 and 1.7 m, and its floor is muddy and covered by mostly eelgrass (Zostera marina) except in the intertidal areas where Manila clam (Ruditapes philippinarum) aquaculture grounds are present [24, 43, 44]. Manila clam and Pacific oyster (Crassostrea gigas) farming are the main aquaculture activities in the lagoon [45]. In Akkeshi Bay, Z. marina occurs in the intertidal to shallow subtidal zones where the depth is less than 2 m, and Z. asiatica dominates in the deeper subtidal zone where the depth is limited to 5 m below the mean low water [46].

A total of 6 study sites were established on the shallow subtidal bottom of Akkeshi-ko Lagoon (BK: at the Bekanbeushi River mouth; CK: Chikarakotan; CL: Central Lagoon; HN: Horonitai; SL: South Lagoon; TB: Toubai), and 2 sites were established in Akkeshi Bay (AK: Aininkappu; SR: Shinryu) (Table 1; Fig 1). See [44] for the detailed explanation of each site and Table 1 for the physical parameters.

Saroma.

The Saroma-ko Lagoon is located on the northeastern part of Hokkaido and is connected to the Sea of Okhotsk by one channel in the east and another channel in the west (Table 1; Fig 1). The area of the lagoon is 152 km2 and has an average depth of 14 m [30]. There are 8 major rivers that flow into the lagoon, with seasonal differences in the freshwater discharges, and the Saromabetsu River in the eastern part of the lagoon is the largest in terms of catchment size [47]. The watersheds of these rivers mainly consist of forests and agricultural lands. Like the Akkeshi-ko Lagoon, the mixture of seawater and freshwater, as well as the presence of water currents within the lagoon, creates a complex physical environment [48], yet the nutrient and salinity gradients within the lagoon are not as notable as the ones in Akkeshi-ko lagoon.

Eelgrass beds are found in the shallow subtidal areas near the coast of the Saroma-ko Lagoon, together with other seagrass species, Z. caspitosa and Z. japonica ([30, 49], personal observation M. Namba). The bottom of the lagoon is mainly sand. Scallop (Mizuhopecten yessoensis) and oyster (Crassostrea gigas) farming are the major primary industries in the lagoon. However, hypoxic events due to the excessive inputs of organic matter from the scallop farming sites and the rivers have been occurring in recent years [48] and are considered to be responsible for the decline of the eelgrass beds over the past decade [30].

A total of 7 study sites (SA1, 2, 3, 5, 6, 7, and 8) were established in the subtidal part of the Saroma-ko Lagoon (Table 1; Fig 1), and these sites corresponded to the periodical monitoring sites by the Hokkaido Aquaculture Promotion Corporation [50].

Notoro.

The Notoro-ko Lagoon is located in the south of Saroma-ko lagoon on the Sea of Okhotsk (Fig 1) and has a maximum depth of 20 m [51] and an area of 58.4 km2 [52]. It is connected to the ocean through a channel that is 324 m side and 13 m deep, and the water exchange rate between the lagoon and the ocean is influenced by tides [51]. Although there are 11 rivers that run into the lagoon, the outfall from Ubaranai River, which is the largest, constitutes only 0.01% of the total water volume of the lagoon due to the small catchment area [52]. Thus, the effects of the rivers are much smaller than that of the ocean, and this makes both salinity and nutrient concentrations almost uniform within the lagoon (personal communication S. Chiba). The salinity of the lagoon is approximately 33, which is similar to the salinity of the Sea of Okhotsk [53]. The watershed of the Notoro-ko Lagoon is smaller than those of the Akkeshi-ko and Saroma-ko Lagoons, and the environment surrounding the lagoon is mainly agricultural lands, with some forests present on the northeastern shore.

In the Notoro-ko Lagoon, two species of seagrasses, Z. marina and Z. caespitosa, form beds on the sandy bottoms of the shallow subtidal areas [21], and Z. japonica occurs in the intertidal areas. The eelgrass beds provide habitat for the shrimp Pandalus latirostris, which is the main fishery resource in the lagoon [21].

A total of 5 study sites were established in the shallow subtidal part of the lagoon (NO1, 2, 3, 4, 5, Table 1; Fig 1). A site was not established in the middle of the lagoon as it is deeper than the depth limit of seagrasses.

Field sampling

In the three lagoons, the first set of field sampling was carried out between June 23 and July 13, 2016 (hereafter ‘summer’) when the eelgrass beds are most productive in this region [11, 54]. The second set of field sampling was undertaken between August 16 and September 28, 2016 (hereafter ‘fall’). Throughout the study, the data obtained from these two seasons were analyzed separately as temporal replicates to assess the reproducibility of the spatial variation. The sampling was conducted during the daytime neap tide so that the effects of tidal currents were minimized. In the summer, SA2 in the Saroma-ko Lagoon was not accessible due to the excess amount of freshwater discharge and subsequent turbidity. In the fall, samples were not collected at SA5 and SA7 in the Saroma-ko Lagoon because the eelgrass at these sites had been decimated by a storm. The water depth in Saroma was deeper than in the other two lagoons, and the depth at SA5 was approximately 3 m. Because of the lack of eelgrass beds in the fall at SA5 where the depth was the deepest, the differences among the average depths in Saroma, Notoro, and Akkeshi were smaller in the fall.

We used memory sensors (AAQ-175 RINKO: JFE Advantech Co. Ltd., Japan, and RINKO-Profiler ASTD102-ALC-R02: JFE Advantech Co., Ltd., Japan) to measure the water temperature, salinity, and water depth at each site. For each site, water samples were collected using a plastic bucket that was washed three times prior to sampling. The water used for the Chl-a measurements was collected using a 138.5 ml darkened polyethylene bottle (one sample per site). For nutrient analysis, 200 ml opaque polyethylene bottles were used to collect water (one sample per site). All the bottles containing water samples were brought back to the laboratory in a darkened cooler box filled with ice.

At each site, we haphazardly collected five replicate samples of mobile invertebrates (hereafter ‘epifauna’) from the aboveground parts of eelgrass using a mesh bag (20 cm diameter, mesh size 0.1 mm) in a circular 0.0314 m2 area in the middle of the eelgrass beds [44]. From these samples, we obtained data on the biomass for both epifauna and eelgrass and then extrapolated the values for 1m2 area from the data. Microalgal samples were obtained by collecting five replicate samples of eelgrass shoots per site, and each sample was placed in a separate plastic zip bag [44] and stored in a darkened cooler box until analysis in the laboratory.

Laboratory procedures

The water sample used for nutrient analysis was filtered through a 0.45 μm nylon membrane filter (Millex—HN Filter Unit, Merck Millipore Ltd., Tullagreen, Carrigtwohill, Ireland) into glass vials and frozen at—20 °C until analysis. The nutrient concentrations, including the amounts of total nitrogen (TN, the sum of NO2-N, NO3-N, and NH4-N) and PO4-P in the water filtrates, were measured simultaneously using an AutoAnalyzer (QuAAtro 39, BL tec, Osaka, Japan).

The non-acidification method of Welschmeyer [55] was used to analyze the Chl-a concentrations. For the water column Chl-a analysis, each water sample was filtered through GF/F glass-fiber filters (Whatman International Ltd., Maidstone, UK). To determine the abundance of epiphytic microalgae from the Chl-a concentration, microalgae on each eelgrass sample was scraped off using a glass slide and filtered through the GF/F filters. The filters were then extracted in 6 ml N,N-dimethylformamide for more than one day and stored at—20 °C in the dark until analysis [56]. The Chl-a concentration within each extract was then measured by using a fluorometer (10-AU-005-CE Fluorometer, Turner Designs, Sunnyvale, CA, USA).

The epifaunal invertebrates collected in the mesh bags were scraped from the eelgrass and filtered through a 0.5 mm sieve. The epifaunal samples retained on the sieve were then fixed with 70% ethanol, and only those with the sizes between 0.5 mm and 8 mm were counted, identified to the lowest possible taxonomic group with a dissecting microscope, and grouped into different size groups for the estimation of the ash-free dry weight using the formulas provided by Edgar [57]. In total, 19 groups of epifauna were identified to the order, class, or phylum level based on the taxonomic knowledge and information available. All eelgrass shoots were then put in small aluminum foil bags and dried at 60 °C for 48 hr or until completely dried in an oven, and they were weighed to obtain the dry mass.

Statistical analysis

The statistical program R [58] was used for all statistical analyses (see S1 Dataset for the all data used for the analyses and S1 Text for the R commands). To assess the spatial variation in the biomasses of the different functional groups in the eelgrass bed community and determine the effects of different spatial scales, a one-factor nested analysis of variance with the sites nested in the lagoons was performed for each season and for each of the following variables; eelgrass (g dry weight per unit area: g DW m-2), microalgae (g Chl-a per unit area: g Chl-a m-2) and epifauna (g ash-free dry weight per unit area: g AFDW m-2). All variables were log-transformed to meet the assumptions of normality and homogeneity. The test was followed by partitioning the variance components using VCAinference code in the VCA package [59].

To determine the abiotic and biotic factors responsible for the variation in the three variables, linear mixed models (LMMs) with Gaussian distributions were used [60]. For the fixed factors, a list of abiotic and biotic candidate variables (water temperature, PO4-P, TN, depth, salinity, microalgae Chl-a, water Chl-a, eelgrass biomass) that could affect the eelgrass biomass, abundance of microalgae and epifaunal biomass was made based on previous observations and the literature (e.g. [24, 44]). Collinearity among the variables was then checked by calculating Pearson’s correlation coefficients for all pairs. Any variables that had absolute values of the coefficient greater than 0.7 [61] were removed from the list, which resulted in a list of five explanatory variables (TN, depth, water temperature, salinity, and water Chl-a) used for the assessment of eelgrass biomass. Microalgae Chl-a and PO4-P were removed from the list due to observed collinearity between salinity and TN, respectively. Different variable combinations were used for the microalgae biomass (TN, depth, water temperature, and water Chl-a, eelgrass biomass) and the epifaunal biomass (TN, depth, water temperature, water Chl-a, eelgrass biomass, microalgae Chl-a) after checking for collinearity among the variables by the abovementioned method. LMM analyses for each of the three functional groups were first carried out to examine all lagoons together (with three random factors: sites nested in the lagoon, lagoon, and season) and then for each lagoon separately (with two random factors: site and season).

The LMMs were fit with the lmer function in the lme4 package [62], and their p-values were obtained by using the lmerTest package [63]. The best model for each of the three producers was chosen based on the AICc, which is Akaike’s information criterion (AIC) adjusted for small sample sizes, where AIC and AICc will be equal at large sample sizes [64]. The AICc was based on restricted maximum likelihood (REML), as the selection was made among the models with nested random factors [65]. To determine the components of the variances explained (R2) by the fixed factors (marginal R2) and the combination of the fixed and random factors (conditional R2) [66], the r.squaredGLMM function in the MuMIn package [67] was used. The variables chosen by the best models were then plotted with linear regression lines shown only if there was a significant relationship between the two variables (linear models, P ≤ 0.05).

Results

Environmental variables

We expected to see the differences in environmental gradients among the lagoons due to water currents at larger spatial scales and the amount of river discharges at finer spatial scales. The results showed that the difference in water temperature was observed among the lagoons, and that variations in salinity and nutrient concentration were most notable within Akkeshi where the amount of river discharge was largest among the three lagoons.

The water temperature in Akkeshi was lower overall than in Saroma and Notoro in the summer, but it was warmer than in the two lagoons on the Okhotsk coast in the fall (Table 1). The average salinity was lower in Akkeshi than in Saroma and Notoro. The salinity varied from 12.3 to 30.9 among the sites in Akkeshi, but it was homogeneous among the sites in Saroma and Notoro (Table 1). Four components of the nutrient concentrations varied among ecosystems and among sites within each ecosystem, and the water nutrient concentrations were higher in Akkeshi than Saroma and Notoro (Fig 2). In addition, the measured concentrations were more heterogeneous among the sites in Akkeshi than in Saroma or Notoro (S1 Table). The water-column Chl-a concentration was 2 to 3-fold higher in Akkeshi than in the two lagoons on the Okhotsk coast during both sampling times (Fig 2), and the concentrations in Saroma and Notoro were relatively homogeneous compared to those in Akkeshi.

thumbnail
Fig 2. Mean (+ SD) concentrations of water nutrients and Chl-a concentrations in the water column.

Nutrients include PO4, NO2, NO3, and NH4. n = 8 for Akkeshi in the summer and fall, n = 6 for Saroma in the summer, n = 5 for Saroma in the fall, n = 5 for Notoro in the summer and fall.

https://doi.org/10.1371/journal.pone.0201791.g002

Eelgrass and microalgae

The above-ground dry weight of eelgrass differed among the lagoons (two-way nested ANOVA: summer F2, 74 = 9.49, P < 0.001; fall F2, 72 = 11.16, P < 0.001) and among the sites within the lagoons (summer F16, 74 = 6.38, P < 0.001; fall F15, 72 = 4.88, P < 0.001) in both the summer and fall (Table 2; Fig 3). The partitioning of the variance components for the nested ANOVA results showed that 50 and 39% of the variances were attributed to the within-lagoon differences, 45 and 50% were attributed to the within-site differences, and only 5 and 10% of the variances were attributed to the among-lagoon differences in the summer and fall, respectively (Table 3).

thumbnail
Fig 3. Mean (+ SD) biomass of eelgrass, microalgae, and epifauna.

n = 8 for Akkeshi in the summer and fall, n = 6 for Saroma in the summer, n = 5 for Saroma in the fall, n = 5 for Notoro in the summer and fall of the above-ground dry weight of eelgrass (g m-2), microalgae weight (g Chl-a m-2), and ash-free dry weight of epifauna (g m-2) at Akkeshi, Saroma, and Notoro for the two temporal replicates (summer and fall). The error bars show the standard deviations.

https://doi.org/10.1371/journal.pone.0201791.g003

thumbnail
Table 2. Mean above-ground dry weight (n = 5) of eelgrass, microalgae, and epifauna.

https://doi.org/10.1371/journal.pone.0201791.t002

thumbnail
Table 3. Variance components (in %) assigned to the different levels by the nested ANOVA.

https://doi.org/10.1371/journal.pone.0201791.t003

The amount of microalgae biomass differed significantly among the lagoons (two-way nested ANOVA: summer F2, 74 = 80.62, P < 0.001; fall F2, 71 = 233.74, P < 0.001), showing that the abundance of microalgae was much higher in Akkeshi (> 100-fold) than in the two Okhotsk lagoons (Table 2; Fig 3). The differences among the sites (two-way nested ANOVA: summer F16, 74 = 7.84, P < 0.001; fall F15, 71 = 12.36, P < 0.001) were also prominent, and the overall variations were higher in Akkeshi than in the other two lagoons. In contrast to the results for eelgrass, the partitioning of the variance components indicated that 50 and 70% of the variances were attributed to the among-lagoon differences, 29 and 21% of the variances were attributed to the within-lagoon differences, and 21 and 9% of the variances were attributed to the within-site differences in the summer and fall, respectively (Table 3).

Epifauna

A total of 104 taxa belonging to 19 taxonomic groups were collected in this study. Gammaridea were most common in the summer in Akkeshi, followed by gastropoda. In the fall, gastropoda was dominant in most sites in Akkeshi. In the summer in Saroma, caprellidea, bivalva, and gastropoda were abundant, whereas gastropoda was the dominant group in the fall. In Notoro, gastropoda was the group that contributed most to the total biomass both in the summer and fall. See S1 Fig for the taxonomic composition of the epifauna for each lagoon and season.

The epifaunal biomass significantly differed among the lagoons (two-way nested ANOVA: summer F2, 76 = 13.79, P < 0.001; fall F2, 72 = 35.02, P < 0.001) for both seasons (Table 2; Fig 3). The lowest epifaunal biomass was observed at the sites of Saroma. There were also differences among the sites within the lagoons (two-way nested ANOVA: summer F16, 76 = 5.00, P < 0.001; fall F15, 72 = 10.66, P < 0.001). BK in Akkeshi had almost 30 times higher epifaunal biomass than the rest of the sites in Akkeshi and all sites in Notoro and Saroma. The partitioning of the variance components showed that 14 and 22% of the variances were from the among-lagoon variations, 38 and 51% were from the within-lagoon variations, and 48 and 27% were from the within-site variations in the summer and fall, respectively (Table 3).

Factors responsible for the spatial variation

For the eelgrass biomass and microalgal abundance, the combinations of explanatory environmental factors differed in the best models (Table 4). Eelgrass was best described by depth (AICc 390.3) when all data from the three lagoons were combined, indicating that the biomass of eelgrass increases with water depth (Fig 4). The fixed factor explained 9% of the total variance of the model, and 54% of the variance was explained when the fixed factor was combined with the random factors (lagoon, site, and season). When the model was analyzed for each lagoon separately, the combinations of the fixed factors differed among lagoons. TN and salinity best described the eelgrass biomass in Akkeshi (AICc 156.5), the water Chl-a concentration and salinity best described the eelgrass biomass in Saroma (AICc 133.7), and depth and salinity best described the eelgrass biomass in Notoro (AICc 59.3).

thumbnail
Fig 4. Scatterplots showing the relationships between eelgrass (above-ground dry weight, g m-2) and the explanatory variables.

The line was obtained from linear regression (only shown if the relationship between the two variables was significant; linear model, P ≤ 0.05), and the values for the explanatory variables (P ≤ 0.05) were chosen in the best LMM for all lagoons combined, Akkeshi, Saroma, and Notoro.

https://doi.org/10.1371/journal.pone.0201791.g004

thumbnail
Table 4. Results of the LMMs for the best combination of the environmental factors responsible for the variation in each producer group (eelgrass, microalgae, and epifauna).

https://doi.org/10.1371/journal.pone.0201791.t004

The microalgae abundance was best described by the combination of water temperature and eelgrass biomass (AICc 468.1), where the microalgae abundance increased with water temperature and eelgrass biomass (Fig 5). The two fixed factors together explained 7% of the total variance of the model, and 81% of the variance was explained by the fixed factors and three random factors. The best models analyzed for each of the three lagoons selected eelgrass biomass and water Chl-a for Akkeshi (AICc 164.1), depth and water Chl-a for Saroma (AICc 144.3), and eelgrass biomass and water temperature for Notoro (AICc 141.4).

thumbnail
Fig 5. Scatterplots showing the relationships between microalgae weight (Chl-a g m-2) and the explanatory variables.

The line was obtained from linear regression (only shown if the relationship between the two variables was significant; linear model, P ≤ 0.05), and the values for the explanatory variables (P ≤ 0.05) were chosen in the best LMM for all lagoons combined, Akkeshi, Saroma, and Notoro.

https://doi.org/10.1371/journal.pone.0201791.g005

For epifauna, only the above-ground dry weight of eelgrass was selected as the explanatory variable of the best model (AICc 413.5; Table 4). This result showed that the epifaunal biomass increased with eelgrass biomass (Fig 6). The above-ground dry weight of eelgrass explained only 4% of the variance, whereas 57% of the total variance was explained by the three random factors and the fixed factor. Although not significant, the microalgae abundance was selected for the best model (AICc 170.4) for Akkeshi. TN and water temperature were selected for Saroma (AICc 150.4), and eelgrass biomass was selected for Notoro (AICc 78.3).

thumbnail
Fig 6. Scatterplots showing the relationships between epifauna (ash-free dry weight, g m-2) and the explanatory variables.

The line was obtained from linear regression (only shown if the relationship between the two variables was significant; linear model, P ≤ 0.05), and the values for the explanatory variables (P ≤ 0.05) were chosen in the best LMM for all lagoons combined, Akkeshi, Saroma, and Notoro.

https://doi.org/10.1371/journal.pone.0201791.g006

Discussion

In this study, we assessed 1) the spatial scale at which eelgrass, epiphytic microalgae and epifauna vary the most, and 2) the abiotic/biotic factors that affect the observed spatial variability in producers and consumers in the eelgrass beds of eastern Hokkaido based on our hypotheses that that susceptibility of each functional group to these large and fine scale gradients is responsible for the spatial scale in which the differences in the biomass are the largest. We found that the amount of variation assigned to different spatial scales varied greatly among the different functional groups. For eelgrass and epifaunal biomass, the within-lagoon variation was greater than the among-lagoon variation. The factors that were responsible for the variation were related to fine scale environmental variations observed within each lagoon. In contrast, the differences in microalgae abundance were more prominent among lagoons than within lagoon. Furthermore, the combinations of responsible factors also varied among functional groups, suggesting that the processes causing the observed variations are different.

Eelgrass

For eelgrass, the variables for explaining the variation in biomass are related to within-lagoon environmental gradients. Depth was selected as one of the predictor variables in the optimal models for each lagoon where eelgrass biomass was greater in deeper depths, which may be related to the taller canopy heights in deeper habitats. For example, SA5 was the deepest site in the Saroma-ko Lagoon, and the biomass was also highest at this site in the lagoon. For eelgrass, increased water depth limits light availability in the deeper parts of the bed, whereas eelgrass is subject to more intense stresses and disturbances in the shallower parts including the intertidal zone [68, 69]. Due to the constraints at both ends of the depth gradient, eelgrass biomass is generally maximized in the intermediate zone of its depth distribution [68, 70]. The observed pattern in this study indicated that the increase in eelgrass biomass occurred from the shallower limits to the optimum depth because our study sites were in the shallower areas of the subtidal zones where it is unlikely that light availability limits plant growth. Nevertheless, depth had overall a weak contribution in the model for all lagoons combined, and this indicates that it is not the most influential factor to explain the observed spatial variations.

Other responsible variables explaining the variation in eelgrass biomass were nutrient concentrations, the amount of phytoplankton in the water column (water Chl-a), and salinity, and these are also related to the environmental gradients within the lagoons. TN was a significant factor only in Akkeshi where nutrient-rich terrestrial water flows from the Bekanbeushi River [71]. In addition to the terrestrial sources, excess nutrients come from shellfish aquaculture [45], and the intensive oyster farming in Akkeshi may somewhat contribute to the presence of the within-lagoon nutrient gradients. The depletion of marine-originated nutrients at the time of eelgrass growth [37, 39] and the lack of both large river inflow and congregated aquacultural activities would explain why TN was not a significant factor in the two Okhotsk lagoons. The amount of water Chl-a is influenced by the water column nutrient concentration and can cause light attenuation in eelgrass bed ecosystems especially under eutrophic conditions [72]. This competition for the light resources between eelgrass and water column phytoplankton could be related to our observed pattern, and we could expand it to understand indirect negative effects of nutrients on eelgrass in relation to increase in water Chl-a for further research on the eelgrass bed ecosystem.

Salinity was selected as a predictor variable that was related to the variation in eelgrass biomass in each of the three lagoons. Similar to the observed nutrient input pattern, the salinity gradient within lagoon was most prominent in Akkeshi due to the presence of the large freshwater input volume from the Bekanbeushi River [24]. In contrast, Saroma and Notoro exhibited less obvious salinity gradients because of the strong and uniform tidal influences across the lagoons. Eelgrass can withstand a wide salinity range [73], but the levels of tolerance vary among populations [74]. More investigations are however required to understand the relationship between eelgrass biomass and salinity gradients.

Microalgae

Unlike eelgrass, the results showed that the difference between the microalgae abundance in Akkeshi and the two sites along the Okhotsk coast (Saroma and Notoro) was great (up to 100-fold higher in Akkeshi), suggesting that the among-lagoon variation was greater than the within-lagoon variation. This difference might be related to the effects of the different ocean currents flowing along the two coasts. The water mass from the Oyashio cold current contained higher nutrient concentrations than that from the Soya warm current during the sampling seasons [36, 37, 39], which was consistent with our observations in the eelgrass beds. Also, the difference in water temperature is created by these two water currents, which is possibly related to the relationship between microalgae and water temperature observed from our results. It is possible that the higher nutrient concentrations along the Oyashio-influenced Pacific coast sustain the higher abundances, or biomass, of microalgae than along the Soya-influenced Okhotsk coast through various mechanisms such as supporting higher turnover rates [75]. Nevertheless, the effect of water temperature on microalgae remains unclear from this study and that the nutrient content was not selected as a predictor in our models, possibly because the nutrients from the ocean were already depleted at the time of sampling, and the rate of nutrient uptake by microalgae is faster than that by eelgrass [76].

Although microalgae biomass was best explained by the among-lagoon variation, there were still some effects of environmental variations observed within lagoons. Eelgrass biomass was one of the factors explaining the microalgal abundance in all lagoons except for Saroma. Aboveground eelgrass tissues provide substrata for the attachment of epiphytic microalgae [77], and eelgrass leaf emergence rate is also related to epiphytic algae load [78]. In addition, the negative relationship between the microalgal abundance and water Chl-a in two of the studied lagoons were observed. As both functional groups have similar light and nutrient requirements [79], the observed pattern may reflect the competitions for these resources among epiphytic microalgae and water column phytoplankton. In Akkeshi, on the other hand, a large portion of the Chl-a in the water comes from benthic microalgae that are detached and resuspended in the water column and not from phytoplankton [43], and this explain the slight difference in the trend.

Epifauna

Our variation partitioning results suggest that scale-specific variation patterns of epifauna are similar to those of eelgrass, in which the spatial variation at the within-lagoon level was much more pronounced than that at the among-lagoon level. The LMM for all lagoons combined, as well as that for Notoro, showed a relationship between epifaunal biomass and eelgrass biomass, suggesting a bottom-up regulation of the abundance of the animals as recorded in previous studies [44, 80]. Gustafsson and Boström [81] reported that the increase in plant biomass that is followed by an increase in seagrass species richness is related to the increase in the associated faunal abundance. Leopardas et al. [82] presented the effect of aboveground seagrass biomass on epifaunal species composition. The ecosystem functions of an eelgrass bed include both habitat provisioning and the supply of food resources to the epifauna [1, 9, 83]. The epifaunal invertebrate grazers studied here consisted mostly of gammaridea, caprellidea, and gastropoda that directly consume epiphytic microalgae on eelgrass blades [84, 85]. For these species, it is likely that the eelgrass biomass provides important habitat space and foraging ground [8688], which could explain the observed pattern regarding to eelgrass and epifaunal biomass and the similarity in the spatial variance pattern over multiple scales. As we see from the interaction between microalgae and eelgrass, and the fact that many epifauna feed on microalgae, there are possible implications for a trophic interaction in the studied eelgrass bed ecosystem. Nevertheless, our result does not show a clear link between microalgae biomass and epifauna biomass, and the observed pattern suggests that the plant-animal interactions in the studied eelgrass beds are more influenced by non-trophic interactions, such as habitat provisioning, than trophic interactions.

There were major differences in the taxonomic composition of epifauna within and among the lagoons, and among the temporal replicates. Although the results of the models for all lagoons combined did not show any effects of abiotic factors on the total epifaunal biomass, we found that TN and water temperature were affecting the epifauna in some degree, which may be explained by the differences in the compositions and functions of the epifaunal communities. For example, nutrient gradients caused by regional-scale eutrophication events [89] as well as regional differences in water temperature [90] are some of the known causes of variations in species composition.

Conclusion

This research fulfills the needs for a local-scale assessment of the functions and services provided by coastal marine ecosystems by examining the spatial variation in the biomasses of different functional groups of plants and animals in eelgrass beds found within three lagoons with unique environmental properties. As expected, differences in the amount of freshwater input and terrestrial runoff from the river input to the lagoons and the current regimes create both among- and within-lagoon environmental variations. Moreover, different patterns of spatial-scale dependency in the variation in the different functional groups, and different combinations of predictors explained the patterns of variation in these organisms were observed from this study. Specifically, we found that those functional groups that exhibit more within- than among-lagoon variations in biomass are influenced by finer scale environmental gradients, and the group that has higher among-lagoon variation is more affected by large scale environmental factors. Thus, it is important to consider and analyze abiotic/biotic variables at multiple spatial scales when assessing the ecosystem functions and services of coastal and estuarine ecosystems, including seagrass beds. Also, future studies should take into account of the differences in environmental properties among various lagoons and other study areas as a same ecosystem can be affected by different combinations of environmental factors based on the properties of the surrounding habitats.

Supporting information

S1 Dataset. The environmental factors and biomass of eelgrass, microalgae, and epifauna.

https://doi.org/10.1371/journal.pone.0201791.s001

(XLSX)

S1 Text. R commands used for the statistical analyses.

https://doi.org/10.1371/journal.pone.0201791.s002

(DOCX)

S1 Table. Water column nutrients (μM) and water Chl-a concentrations (μg/l) in the summer and fall 2016 at the study sites in Akkeshi, Saroma, and Notoro.

https://doi.org/10.1371/journal.pone.0201791.s003

(DOCX)

S1 Fig. Mean biomass of epifauna (g m-2; n = 5).

The biomass is expressed by 9 major taxa for each site in Akkeshi, Saroma, and Notoro in the summer and fall. Minor taxonomical groups were grouped to one category (shown as Others). The symbols (*) on the graphs indicate no data.

https://doi.org/10.1371/journal.pone.0201791.s004

(PDF)

Acknowledgments

We would like to thank T. Isada, K. Miyashita, T. Noda, H. Glon, K. Sudo and three anonymous reviewers for improving the manuscript and S. Hamano and H. Katsuragawa for their technical support. We are also thankful for the assistance provided by T. Hasegawa, M. Hashimoto, D. Hayashi, N. Hookabe, T. Imamura, M. Ito, M. Kishimoto, K. Momota, C. Smith, K. Takahoshi, Z. Tamura, N. Ueno, R. Yamamoto, and T. Yorisue during the field and laboratory procedures. K. Sakaguchi and other members of the Aquaculture Fishery Cooperative of Saroma Lake, T. Kawajiri and other members of the Nishi-Abashiri Fisheries Union, and the Abashiri City Fisheries Science Center and Abashiri City provided logistical support during our field work in Saroma and Notoro. This research was partially funded by the Environment Research and Technology Development Fund (S-15 Predicting and Assessing Natural Capital and Ecosystem Services (PANCES)) of the Ministry of the Environment, Japan, and a grant from the Aquaculture Fishery Cooperative of Saroma Lake.

References

  1. 1. Barbier EB, Hacker SD, Kennedy C, Koch EW, Stier AC, Silliman BR. The value of estuarine and coastal ecosystem services. Ecol Monogr. 2010;81: 169–193.
  2. 2. Seitz RD, Wennhage H, Bergström U, Lipcius RN, Ysebaert T. Ecological value of coastal habitats for commercially and ecologically important species. ICES J Mar Sci J du Cons. 2014;71: 648–665.
  3. 3. Costanza R, D’Arge R, de Groot R, Farber S, Grasso M, Hannon B, et al. The value of the world’s ecosystem services and natural capital. Nature. 1997;387: 253–260.
  4. 4. Short FT, Neckles HA. The effects of global climate change on seagrasses. Aquat Bot. 1999;63: 169–196.
  5. 5. de Groot RS, Wilson MA, Boumans RMJ. A typology for the classification, description and valuation of ecosystem functions, goods and services. Ecol Econ. 2002;41: 393–408.
  6. 6. Duffy JE. Biodiversity and the functioning of seagrass ecosystems. Mar Ecol Prog Ser. 2006;311: 233–250.
  7. 7. Duarte CM, Chiscano CL. Seagrass biomass and production: a reassessment. Aquat Bot. 1999;65: 159–174. Available: http://www.sciencedirect.com/science/article/pii/S0304377099000388
  8. 8. Moncreiff CA, Sullivan MJ. Trophic importance of epiphytic algae in subtropical seagrass beds: evidence from multiple stable isotope analyses. Mar Ecol Prog Ser. 2001;215: 93–106.
  9. 9. Nakaoka M. Plant–animal interactions in seagrass beds: ongoing and future challenges for understanding population and community dynamics. Popul Ecol. 2005;47: 167–177.
  10. 10. Burkholder JM, Tomasko DA, Touchette BW. Seagrasses and eutrophication. J Exp Mar Bio Ecol. 2007;350: 46–72.
  11. 11. Hasegawa N, Hori M, Mukai H. Seasonal shifts in seagrass bed primary producers in a cold-temperate estuary: Dynamics of eelgrass Zostera marina and associated epiphytic algae. Aquat Bot. 2007;86: 337–345.
  12. 12. Kollars NM, Henry AK, Whalen MA, Boyer KE, Cusson M, Eklöf JS, et al. Meta-Analysis of Reciprocal Linkages between Temperate Seagrasses and Waterfowl with Implications for Conservation. Frontiers in plant science. 2017;8. pmid:29312384
  13. 13. Fredette TJ, Diaz RJ, van Montfrans J, Orth RJ. Secondary production within a seagrass bed (Zostera marina and Ruppia maritima) in lower Chesapeake Bay. Estuaries. 1990;13: 431.
  14. 14. Duffy JE, Richardson JP, France KE. Ecosystem consequences of diversity depend on food chain length in estuarine vegetation. Ecol Lett. 2005;8: 301–309.
  15. 15. Hemminga MA, Duarte CM. Seagrass Ecology. Cambridge, United Kingdom: Cambridge University Press; 2000.
  16. 16. Duffy JE, Hughes AR, Moksnes P-O. Ecology of seagrass communities. In: Bertness MD, Bruno JF, Silliman BR, Stachowicz JJ, editors. Marine community ecology and conservation. Sunderland, MA: Sinauer; 2013. pp. 271–297.
  17. 17. Davis BC, Fourqurean JW. Competition between the tropical alga, Halimeda incrassata, and the seagrass, Thalassia testudinum. Aquat Bot. 2001;71: 217–232.
  18. 18. Duffy JE, Reynolds PL, Boström C, Coyer JA, Cusson M, Donadi S, et al. Biodiversity mediates top–down control in eelgrass ecosystems: a global comparative-experimental approach. Ecol Lett. 2015;18: 696–705. Available: http://onlinelibrary.wiley.com/doi/10.1111/ele.12448/full pmid:25983129
  19. 19. Boström C, Bonsdorff E. Zoobenthic community establishment and habitat complexity-the importance of seagrass shoot-density, morphology and physical disturbance for faunal recruitment. Mar Ecol Prog Ser. 2000;205: 123–138. Available: https://www.abo.fi/personal/en/media/16577/bostrombonsdorffmeps2000.pdf
  20. 20. Boström C, Jackson EL, Simenstad CA. Seagrass landscapes and their effects on associated fauna: A review. Estuar Coast Shelf Sci. 2006;68: 383–403.
  21. 21. Shinomiya Y, Chiba S, Kanamori M, Hashizume S, Yoshino K, Goshima S. Importance of patch size variation for the population persistence of a decapod crustacean in seagrass beds. Mar Ecol Prog Ser. 2017;570: 157–171. Available: http://www.int-res.com/abstracts/meps/v570/p157-171/
  22. 22. Orth RJ, Heck KLJ, van Montfrans J. Faunal communities in seagrass beds: A review of the influence of plant structure and prey characteristics on predator-prey relationships. Estuaries. 1984;7: 339–350.
  23. 23. Ferdie M, Fourqurean JW. Responses of seagrass communities to fertilization along a gradient of relative availability of nitrogen and phosphorus in a carbonate environment. Limnol Oceanogr. 2004;49: 2082–2094.
  24. 24. Yamada K, Hori M, Tanaka Y, Hasegawa N, Nakaoka M. Temporal and spatial macrofaunal community changes along a salinity gradient in seagrass meadows of Akkeshi-ko estuary and Akkeshi Bay, northern Japan. Hydrobiologia. 2007;592: 345–358. http://dx.doi.org/10.1007/s10750-007-0767-6
  25. 25. Schils T, Coppejans E. Spatial variation in subtidal plant communities around the Socotra Archipelago and their biogeographic affinities within the Indian Ocean. Mar Ecol Prog Ser. 2003;251: 103–114.
  26. 26. Blanchette CA, Broitman BR, Gaines SD. Intertidal community structure and oceanographic patterns around Santa Cruz Island, CA, USA. Mar Biol. 2006;149: 689–701. Available: http://link.springer.com/article/10.1007/s00227-005-0239-3
  27. 27. Damme S Van, Struyf E, Maris T, Ysebaert T, Dehairs F, Tackx M, et al. Spatial and temporal patterns of water quality along the estuarine salinity gradient of the Scheldt estuary (Belgium and The Netherlands): results of an integrated monitoring approach. Hydrobiologia. 2005;540: 29–45.
  28. 28. McIver R, Milewski I, Lotze HK. Land use and nitrogen loading in seven estuaries along the southern Gulf of St. Lawrence, Canada. Estuar Coast Shelf Sci. 2015;165: 137–148.
  29. 29. Krause-Jensen D, Pedersen MF, Jensen C. Regulation of eelgrass (Zostera marina) cover along depth gradients in Danish coastal waters. Estuaries. 2003;26: 866–877.
  30. 30. Katsuki K, Seto K, Nomura R, Maekawa K, Khim B-K. Effect of human activity on Lake Saroma (Japan) during the past 150 years: Evidence by variation of diatom assemblages. Estuar Coast Shelf Sci. 2009;81: 215–224.
  31. 31. Palmer TA, Montagna PA, Pollack JB, Kalke RD, Deyoe HR. The role of freshwater inflow in lagoons, rivers, and bays. Hydrobiologia. 2011;667: 49–67. http://dx.doi.org/10.1007/s10750-011-0637-0
  32. 32. Maranon E, Holligan PM, Varela M, Mourino B, Bale AJ. Basin-scale variability of phytoplankton biomass, production and growth in the Atlantic Ocean. Deep Res Part I-Oceanographic Res Pap. 2000;47: 825–857.
  33. 33. Mizushima T. Ecology and catch of the prawn, Pandalus kessleri in Zostera beds in Notsuke Bay, Hokkaido, Japan. Seaweed Beds. 1981; 57–74.
  34. 34. Seto M, Sakurai I, Matsuoka G, Yamashita T. Habitable wave condition for Pandalus kesseri. Proc Civ Eng Ocean. 1996;12: 463–466.
  35. 35. Biodiversity Center of Japan. The Report for the 7th Natural Environmental Survey Shallow Marine Ecosystem Survey (Aquatic Vegetation Survey) [Internet]. Yamanashi, Japan; 2008. Available: http://www.biodic.go.jp/reports2/6th/6_moba19/6_moba19.pdf
  36. 36. Saito H, Tsuda A, Kasai H. Nutrient and plankton dynamics in the Oyashio region of the western subarctic Pacific Ocean. Deep Sea Res Part II Top Stud Oceanogr. 2002;49: 5463–5486. Available: http://www.sciencedirect.com/science/article/pii/S0967064502002047
  37. 37. Kasai H, Saito H, Yoshimori A, Taguchi S. Variability in timing and magnitude of spring bloom in the Oyashio region, the western subarctic Pacific off Hokkaido, Japan. Fish Oceanogr. 1997;6: 118–129. Available: http://onlinelibrary.wiley.com/doi/10.1046/j.1365-2419.1997.00034.x/full
  38. 38. Kasai H, Nakano Y, Ono T, Tsuda A. Seasonal change of oceanographic conditions and chlorophyll a vertical distribution in the southwestern Okhotsk Sea during the non-iced season. J Oceanogr. 2010;66: 13–26.
  39. 39. Kasai H, Nagata T, Murai K, Katakura S, Tateyama /k, Hamaoka S. Seasonal change in oceanographic environments and the influence of interannual variation in the timing of sea-ice retreat on chlorophyll a concentration in the coastal water of northeastern Hokkaido along the Okhotsk Sea (in Japanese). Bull Coast Oceanogr. 2017;54: 181–192. Available: http://ci.nii.ac.jp/naid/40021120147/
  40. 40. Takizawa T. Characteristics of the Sôya Warm Current in the Okhotsk Sea. J Oceanogr Soc Japan. 1982;38: 281–292.
  41. 41. Cavalieri DJ, Parkinson CL. On the relationship between atmospheric circulation and the fluctuations in the sea ice extents of the Bering and Okhotsk Seas. J Geophys Res Ocean. 1987;92: 7141–7162.
  42. 42. Iizumi H, Taguchi S, Minami T, Maekawa S. Distribution and variability of nutrients, chlorophyll a, particulate organic matters, and their carbon and nitrogen contents, in Akkeshi-ko, an estuary in northern Japan. Bull Hokkaido Natl Fish Res Inst. 1995;59: 43–67. Available: http://agriknowledge.affrc.go.jp/RN/2010520965.pdf
  43. 43. Kasim M, Mukai H. Contribution of benthic and epiphytic diatoms to clam and oyster production in the Akkeshi-ko estuary. J Oceanogr. 2006;62: 267–281. https://doi.org/10.1007/s10872-006-0051-9
  44. 44. Momota K, Nakaoka M. Influence of different types of sessile epibionts on the community structure of mobile invertebrates in an eelgrass bed. PeerJ. 2017;5: e2592. Available: https://peerj.com/preprints/2591
  45. 45. Abe H. Studies on the role of primary producer diversity and suspension-feeding bivalves in the lower trophic ecosystem in subarctic estuaries and lagoons -integration between field observation and numerical modelling- [dissertation]. Hokkaido University; 2016.
  46. 46. Watanabe M, Nakaoka M, Mukai H. Seasonal variation in vegetative growth and production of the endemic Japanese seagrass Zostera asiatica: a comparison with sympatric Zostera marina. Bot Mar. 2005;48: 266–273.
  47. 47. Sato T, Imazu Y, Sakawa T, Kazama T, Wakabayashi T, Uhrenholdt T. Modeling of integrated marine ecosystem including the generation-tracing type scallop growth model. Ecol Modell. 2007;208: 263–285.
  48. 48. Terasaki E, Morita K, Yasuda M, Maekawa K, Montani S. Spatial distribution of phytopigments and organic matter in surface sediments in Lake Saroma (Hokkaido, Japan). La Mer. 2013;51: 119–128. Available: http://eprints.lib.hokudai.ac.jp/dspace/handle/2115/56666
  49. 49. Onishi Y, Chiba S, Goshima S. Reproductive cycle of the protandric shrimp Pandalus latirostris in Saroma lagoon. Benthos Res. 2001;56: 9–20.
  50. 50. Hokkaido Aquaculture Promotion Corporation. 2015 Saroma-ko lagoon Eelgrass Distribution Survey Report. Sapporo; 2015.
  51. 51. Nishino Y, Satoh T, Shiomoto A, Taniguchi A. Lower trophic levels of coastal marine ecosystem in lagoon Notoroko on the coastal Okhotsk Sea (in Japanese). J Agric Sci Tokyo Univ Agric. 2014;59: 11–20
  52. 52. Nishino Y, Satoh T, Taniguchi A. Water condition in lagoon Notoro-ko in the Eastern Hokkaido −Seasonal changes of temperature, salinity, dissolved oxygen and nutrients−. Eco-Engineering. 2014;26: 3–9. Available: https://www.jstage.jst.go.jp/article/seitaikogaku/26/1/26_3/_pdf
  53. 53. Imada K, Sakazaki S, Kawajiri T, Kobayashi K. Lake basin and saline environments profiles of four brackish lakes in Abashiri, northern Hokkaido (in Japanese). Sci Reports Hokkaido Fish Hatch. 1995; 37–48. Available: http://ci.nii.ac.jp/naid/40004700088
  54. 54. Nakaoka M, Aioi K. Ecology of seagrasses Zostera spp.(Zosteraceae) in Japanese waters: a review. Otsuchi Mar Sci. 2001;26: 7–22. Available: http://repository.dl.itc.u-tokyo.ac.jp/dspace/handle/2261/5188
  55. 55. Welschmeyer NA. Fluorometric analysis of chlorophyll a in the presence of chlorophyll b and pheopigments. Limnol Oceanogr. 1994;39: 1985–1992. Available: http://www.jstor.org/stable/2838404
  56. 56. Suzuki R, Ishimaru T. An improved method for the determination of phytoplankton chlorophyll using N, N-dimethylformamide. J Oceanogr Soc Japan. 1990;46: 190–194. Available: http://www.springerlink.com/index/RN8724Q88T57T067.pdf
  57. 57. Edgar GJ. The use of the size structure of benthic macrofaunal communities to estimate faunal biomass and secondary production. J Exp Mar Bio Ecol. 1990;137: 195–214.
  58. 58. R Core Team. R: A language and environment for statistical computing [Internet]. Vienna; 2017. Available: http://www.r-project.org/
  59. 59. Schuetzenmeister A. Package “VCA” [Internet]. 2016. Available: https://cran.r-project.org/web/packages/VCA/VCA.pdf
  60. 60. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, et al. Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2009;24: 127–135. pmid:19185386
  61. 61. Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography. 2013;36: 27–46. Available: http://onlinelibrary.wiley.com/doi/10.1111/j.1600-0587.2012.07348.x/full
  62. 62. Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, et al. Package “lme4” [Internet]. 2016. Available: http://kambing.ui.ac.id/cran/web/packages/lme4/lme4.pdf
  63. 63. Kuznetsova A, Brockhoff PB, Christensen RHB. Package “lmerTest” [Internet]. 2015. Available: http://cran.uib.no/web/packages/lmerTest/lmerTest.pdf
  64. 64. Burnham KP, Anderson DR. Multimodel inference: understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33: 261–304.
  65. 65. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology. Spring Science and Business Media. New York, NY: Springer Science+Business Media; 2009.
  66. 66. Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4: 133–142. Available: http://onlinelibrary.wiley.com/doi/10.1111/j.2041-210x.2012.00261.x/full
  67. 67. Bartoń K. MuMIn: Multi-model inference. R package version 1.9. 13. Compr R Arch Netw (CRAN), Vienna, Austria. 2013;
  68. 68. Krause-Jensen D, Middleboe AL, Sand-Jensen K, Christensen PB. Eelgrass, Zostera marina, growth along depth gradients: upper boundaries of the variation as a powerful predictive tool. Oikos. 2000;91: 233. Available: http://ezproxy.library.dal.ca/login?url = http://search.ebscohost.com/login.aspx?direct=true&db=eih&AN=8226997&site=ehost-live
  69. 69. Riis T, Hawes I. Effect of wave exposure on vegetation abundance, richness and depth distribution of shallow water plants in a New Zealand lake. Freshw Biol. 2003;48: 75–87.
  70. 70. Chambers PA, Kaiff J. Depth distribution and biomass of submersed aquatic macrophyte communities in relation to Secchi depth. Can J Fish Aquat Sci. 1985;42: 701–709. Available: http://www.nrcresearchpress.com/doi/abs/10.1139/f85-090
  71. 71. Akabane H, Kishi M, Mukai H, Iizumi H. The responses of the ecosystem of Akkeshi Lake, an estuary in Northern Japan, to nutrients input from terrestrial area. Bull Coast Oceanogr. 2003;40: 171–179. Available: http://ci.nii.ac.jp/naid/110007999287/
  72. 72. Olesen B. Regulation of light attenuation and eelgrass Zostera marina depth distribution in a Danish embayment. Marine Ecology Progress Series. 1996;134(1–3):187–94.
  73. 73. Touchette BW. Seagrass-salinity interactions: Physiological mechanisms used by submersed marine angiosperms for a life at sea. J Exp Mar Bio Ecol. 2007;350: 194–215.
  74. 74. Salo T, Pedersen MF, Boström C. Population specific salinity tolerance in eelgrass (Zostera marina). J Exp Mar Bio Ecol. 2014;461: 425–429.
  75. 75. Huete-Ortega M, Rodríguez-Ramos T, López-Sandoval DC, Cermeño P, Blanco JM, Palomino RL, et al. Distinct patterns in the size-scaling of abundance and metabolism in coastal and open-ocean phytoplankton communities. Mar Ecol Prog Ser. 2014;515: 61–71.
  76. 76. Pedersen MF, Borum J. Nutrient control of algal growth in estuarine waters. Nutrient limitation and the importance of nitrogen requirements and nitrogen storage among phytoplankton and species of macroalgae. Mar Ecol Prog Ser. 1996;142: 261–272. Available: http://www.int-res.com/abstracts/meps/v142/p261-272/
  77. 77. Gordon N, Adams JB, Bate GC. Epiphytes of the St. Lucia Estuary and their response to water level and salinity changes during a severe drought. Aquat Bot. 2008;88: 66–76.
  78. 78. Ruesink JL. Epiphyte load and seagrass performance are decoupled in an estuary with low eutrophication risk. J Exp Mar Bio Ecol. 2016;481: 1–8.
  79. 79. Sand-Jensen K, Borum J. Interactions among phytoplankton, periphyton, and macrophytes in temperate fresh-waters and estuaries. Aquat Bot. 1991;41: 137–175.
  80. 80. Edgar GJ. The influence of plant structure on the species richness, biomass and secondary production of macrofaunal assemblages associated with Western Australian seagrass beds. J Exp Mar Bio Ecol. 1990;137: 215–240. Available: http://www.sciencedirect.com/science/article/pii/002209819090186G
  81. 81. Gustafsson C, Boström C. Biodiversity influences ecosystem functioning in aquatic angiosperm communities. Oikos. 2011;120: 1037–1046.
  82. 82. Leopardas V, Uy W, Nakaoka M. Benthic macrofaunal assemblages in multispecific seagrass meadows of the southern Philippines: Variation among vegetation dominated by different seagrass species. J Exp Mar Bio Ecol. 2014;457: 71–80.
  83. 83. Nordlund LM, Koch EW, Barbier EB, Creed JC. Correction: Seagrass ecosystem services and their variability across genera and geographical regions. PLoS One. 2017;12: e0169942. pmid:28056075
  84. 84. Valentine JF, Heck KL Jr. Seagrass herbivory: evidence for the continued grazing of marine grasses. Mar Ecol Prog Ser. 1999;176: 291–302.
  85. 85. Valentine JF, Duffy JE. The central role of grazing in seagrass ecology. In: Larkum AWD, Orth RJ, Duarte CM, editors. Seagrasses: Biology, ecology, and conservation. Dordrecht, The Netherlands: Springer; 2006. pp. 463–501.
  86. 86. Stoner AW. The role of seagrass biomass in the organization of benthic macrofaunal assemblages. Bull Mar Sci. 1980;30: 537–551. Available: http://www.ingentaconnect.com.ezproxy.library.dal.ca/search/download?pub=infobike://umrsmas/bullmar/1980/00000030/00000003/art00001&mimetype=application/pdf
  87. 87. Doi H, Matsumasa M, Fujikawa M, Kanou K, Suzuki T, Kikuchi E. Macroalgae and seagrass contribution to gastropods in sub-tropical and temperate tidal flats. J Mar Biol Assoc United Kingdom. 2009;89: 399–404.
  88. 88. Sturaro N, Lepoint G, Vermeulen S, Gobert S. Multiscale variability of amphipod assemblages in Posidonia oceanica meadows. J Sea Res. 2015;95: 258–271.
  89. 89. Schmidt AL, Coll M, Lotze HK. Regional-scale differences in eutrophication effects on eelgrass-associated (Zostera marina) macrofauna. Estuaries and Coasts. 2017;40: 1096–1112.
  90. 90. Namba M, Lotze HK, Schmidt AL. Large-scale differences in community structure and ecosystem services of eelgrass (Zostera marina) beds across three regions in Eastern Canada. Estuaries and Coasts. 2018;41: 177–192. https://doi.org/10.1007/s12237-017-0271-9